14 results
Transport by deep convection in basin-scale geostrophic circulation: turbulence-resolving simulations
- Catherine A. Vreugdenhil, Bishakhdatta Gayen, Ross W. Griffiths
-
- Journal:
- Journal of Fluid Mechanics / Volume 865 / 25 April 2019
- Published online by Cambridge University Press:
- 26 February 2019, pp. 681-719
-
- Article
- Export citation
-
Direct numerical simulations are used to investigate the nature of fully resolved small-scale convection and its role in large-scale circulation in a rotating $f$-plane rectangular basin with imposed surface temperature difference. The large-scale circulation has a horizontal geostrophic component and a deep vertical overturning. This paper focuses on convective circulation with no wind stress, and buoyancy forcing sufficiently strong to ensure turbulent convection within the thermal boundary layer (horizontal Rayleigh numbers $Ra\approx 10^{12}{-}10^{13}$). The dynamics are found to depend on the value of a convective Rossby number, $Ro_{\unicode[STIX]{x0394}T}$, which represents the strength of buoyancy forcing relative to Coriolis forces. Vertical convection shifts from a mean endwall plume under weak rotation ($Ro_{\unicode[STIX]{x0394}T}>10^{-1}$) to ‘open ocean’ chimney convection plus mean vertical plumes at the side boundaries under strong rotation ($Ro_{\unicode[STIX]{x0394}T}<10^{-1}$). The overall heat throughput, horizontal gyre transport and zonally integrated overturning transport are then consistent with scaling predictions for flow constrained by thermal wind balance in the thermal boundary layer coupled to vertical advection–diffusion balance in the boundary layer. For small Rossby numbers relevant to circulation in an ocean basin, vertical heat transport from the surface layer into the deep interior occurs mostly in ‘open ocean’ chimney convection while most vertical mass transport is against the side boundaries. Both heat throughput and the mean circulation (in geostrophic gyres, boundary currents and overturning) are reduced by geostrophic constraints.
Ablation of sloping ice faces into polar seawater
- Mainak Mondal, Bishakhdatta Gayen, Ross W. Griffiths, Ross C. Kerr
-
- Journal:
- Journal of Fluid Mechanics / Volume 863 / 25 March 2019
- Published online by Cambridge University Press:
- 28 January 2019, pp. 545-571
-
- Article
- Export citation
-
The effects of the slope of an ice–seawater interface on the mechanisms and rate of ablation of the ice by natural convection are examined using turbulence-resolving simulations. Solutions are obtained for ice slopes $\unicode[STIX]{x1D703}=2^{\circ }{-}90^{\circ }$, at a fixed ambient salinity and temperature, chosen to represent common Antarctic ocean conditions. For laminar boundary layers the ablation rate decreases with height, whereas in the turbulent regime the ablation rate is found to be height independent. The simulated laminar ablation rates scale with $(\sin \unicode[STIX]{x1D703})^{1/4}$, whereas in the turbulent regime it follows a $(\sin \unicode[STIX]{x1D703})^{2/3}$ scaling, both consistent with the theoretical predictions developed here. The reduction in the ablation rate with shallower slopes arises as a result of the development of stable density stratification beneath the ice face, which reduces turbulent buoyancy fluxes to the ice. The turbulent kinetic energy budget of the flow shows that, for very steep slopes, both buoyancy and shear production are drivers of turbulence, whereas for shallower slopes shear production becomes the dominant mechanism for sustaining turbulence in the convective boundary layer.
Turbulent horizontal convection under spatially periodic forcing: a regime governed by interior inertia
- Madelaine G. Rosevear, Bishakhdatta Gayen, Ross W. Griffiths
-
- Journal:
- Journal of Fluid Mechanics / Volume 831 / 25 November 2017
- Published online by Cambridge University Press:
- 13 October 2017, pp. 491-523
-
- Article
- Export citation
-
Differential heating applied at a single horizontal boundary forces ‘horizontal convection’, even when there is no net heat flux through the boundary. However, almost all studies of horizontal convection have been limited to a special class of problem in which temperature or heat flux differences were applied in only one direction and over the horizontal length of a box (the Rossby problem; Rossby, Deep-Sea Res., vol. 12, 1965, pp. 9–16). These conditions strongly constrain the flow. Here we report laboratory experiments and direct numerical simulations (DNS) extending the results of Griffiths & Gayen (Phys. Rev. Lett., vol. 115, 2015, 204301) for horizontal convection forced by boundary conditions imposed in a two-dimensional periodic array at a horizontal boundary. The experiments use saline and freshwater fluxes at a permeable base with the imposed boundary salinity having a horizontal length scale one quarter of the width of the box. The flow reaches a state in which the net boundary buoyancy flux vanishes and the bulk of the fluid shows an inertial range of turbulence length scales. A regime transition is seen for increasing water depth, from an array of individual coherent plumes on the forcing scale to convection dominated by emergent larger scales of overturning. The DNS explore the analogous thermally forced case with sinusoidal boundary temperature of wavenumber $n=4$, and are used to examine the Rayleigh number ($Ra$) dependence for shallow- and deep-water cases. For shallow water the flow transitions with increasing $Ra$ from laminar to turbulent boundary layer regimes that are familiar from the Rossby problem and which have normalised heat transport scaling as $Nu\sim Ra^{1/5}$ and $Nu\sim (Ra\,Pr)^{1/5}$, with $Nu$ the Nusselt number and $Pr$ the Prandtl number, in this case maintaining a stable array of coherent turbulent plumes. For deep-water and large $Ra$ the laminar scaling transitions to $Nu\sim (Ra\,Pr)^{1/4}$, with the scales of turbulence extending to the dimensions of the box. The $1/4$ power law regime is explained in terms of the momentum of symmetric, inviscid large scales of motion in the interior coupled to diffusive loss of heat through stabilised parts of the boundary layer. The turbulence production is predominantly by shear instability rather than convection, with viscous dissipation distributed throughout the bulk of the fluid. These conditions are not seen in the highly asymmetric flow in the Rossby problem even at Rayleigh numbers up to six orders of magnitude greater than the transition found here. The new inertial interior regime has the rate of supply of available potential energy, and its removal by mixing of density, increasing as $Ra^{5/4}$, which is faster than $Ra^{6/5}$ in the Rossby problem. Irreversible mixing is confined close to the forcing boundary and is very much larger than the viscous dissipation, which is proportional to $Ra$.
Geostrophic and chimney regimes in rotating horizontal convection with imposed heat flux
- Catherine A. Vreugdenhil, Ross W. Griffiths, Bishakhdatta Gayen
-
- Journal:
- Journal of Fluid Mechanics / Volume 823 / 25 July 2017
- Published online by Cambridge University Press:
- 15 June 2017, pp. 57-99
-
- Article
- Export citation
-
Convection in a rotating rectangular basin with differential thermal forcing at one horizontal boundary is examined using laboratory experiments. The experiments have an imposed heat flux boundary condition, are at large values of the flux Rayleigh number ($Ra_{F}\sim O(10^{13}{-}10^{14})$ based on the box length $L$), use water with Prandtl number $Pr\approx 4$ and have a small depth to length aspect ratio. The results show the conditions for transition from non-rotating horizontal convection governed by an inertial–buoyancy balance in the thermal boundary layer, to circulation governed by geostrophic flow in the boundary layer. The geostrophic balance constrains mean flow and reduces the heat transport as Nusselt number $Nu\sim (Ra_{F}Ro)^{1/6}$, where $Ro=B^{1/2}/f^{3/2}L$ is the convective Rossby number, $B$ is the imposed buoyancy flux and $f$ is the Coriolis parameter. Thus flow in the geostrophic boundary layer regime is governed by the relative roles of horizontal convective accelerations and Coriolis accelerations, or buoyancy and rotation, in the boundary layer. Experimental evidence suggests that for more rapid rotation there is another transition to a regime in which the momentum budget is dominated by fluctuating vertical accelerations in a region of vortical plumes, which we refer to as a ‘chimney’ following related discussion of regions of deep convection in the ocean. Coupling of the chimney convection in the region of destabilising boundary flux to the diffusive boundary layer of horizontal convection in the region of stabilising boundary flux gives heat transport independent of rotation in this ‘inertial chimney’ regime, and the new scaling $Nu\sim Ra_{F}^{1/4}$. Scaling analysis predicts the transition conditions observed in the experiments, as well as a further ‘geostrophic chimney’ regime in which the vertical plumes are controlled by local geostrophy. When $Ro<10^{-1}$, the convection is also observed to produce a set of large basin-scale gyres at all depths in the time-averaged flow.
Simulation of convection at a vertical ice face dissolving into saline water
- Bishakhdatta Gayen, Ross W. Griffiths, Ross C. Kerr
-
- Journal:
- Journal of Fluid Mechanics / Volume 798 / 10 July 2016
- Published online by Cambridge University Press:
- 31 May 2016, pp. 284-298
-
- Article
- Export citation
-
We investigate the convection and dissolution rate generated when a wall of ice dissolves into seawater under Antarctic Ocean conditions. In direct numerical simulations three coupled interface equations are used to solve for interface temperature, salinity and ablation velocity, along with the boundary layer flow and transport. The main focus is on ambient water temperatures between $-1\,^{\circ }\text{C}$ and $6\,^{\circ }\text{C}$ and salinities around 35 ‰, where diffusion of salt to the ice–water interface depresses the freezing point and enhances heat diffusion to the ice. We show that fluxes of both heat and salt to the interface are significant in governing the dissolution of ice, and the ablation velocity agrees well with experiments and a recent theoretical prediction. The same turbulent flow dynamics and ablation rate are expected to apply at any depth in a deeper ocean water column (after choosing the relevant pressure coefficient for the liquidus temperature). At Grashof numbers currently accessible by direct numerical simulation, turbulence is generated both directly from buoyancy flux and from shear production in the buoyancy-driven boundary layer flow, whereas shear production by the convective flow is expected to be more important at geophysical scales. The momentum balance in the boundary layer is dominated by buoyancy forcing and wall stress, with the latter characterised by a large drag coefficient.
Stability transitions and turbulence in horizontal convection
- Bishakhdatta Gayen, Ross W. Griffiths, Graham O. Hughes
-
- Journal:
- Journal of Fluid Mechanics / Volume 751 / 25 July 2014
- Published online by Cambridge University Press:
- 25 June 2014, pp. 698-724
-
- Article
- Export citation
-
Recent results have shown that convection forced by a temperature gradient along one horizontal boundary of a rectangular domain at a large Rayleigh number can be turbulent in parts of the flow field. However, the conditions for onset of turbulence, the dependence of flow and heat transport on Rayleigh number, and the roles of large and small scales in the flow, have not been established. We use three-dimensional direct numerical simulation (DNS) and large-eddy simulation (LES) over a wide range of Rayleigh numbers, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}Ra\sim 10^8\mbox{--}10^{15}$, for Prandtl number $Pr=5$ and a small aspect ratio, and show that a sequence of several stability transitions at $Ra \sim 10^{10}\mbox{--} 10^{11}$ defines a change from laminar to turbulent flow. The Prandtl number dependence too is examined at $Ra = 5.86 \times 10^{11}$. At the smallest $Ra$ considered the thermal boundary layer is characterized by a balance of viscous stress and buoyancy, whereas inertia and buoyancy dominate in the large-$Ra$ regime. The change in the momentum balance is accompanied by turbulent enhancement of the overall heat transfer, although both laminar and turbulent regimes give $Nu\sim Ra^{1/5}$. The results support both viscous and inviscid theoretical scaling models from previous work. The mechanical energy budget for an intermediate range of Rayleigh numbers above onset of instability ($10^{10}<Ra<10^{13}$) reveals that the small scales of motion are produced predominantly by thermal convection, whereas at $Ra \ge 10^{14}$ shear instability of the large-scale flow begins to play a dominant role in sustaining the small-scale turbulence. Extrapolation to ocean conditions requires knowledge of the inertial regime identified here, but the simulations show that the corresponding asymptotic balance has not been fully realized by $Ra \sim 10^{15}$.
Available potential energy in Rayleigh–Bénard convection
- Graham O. Hughes, Bishakhdatta Gayen, Ross W. Griffiths
-
- Journal:
- Journal of Fluid Mechanics / Volume 729 / 25 August 2013
- Published online by Cambridge University Press:
- 09 August 2013, R3
-
- Article
- Export citation
-
The mechanical energy budget for thermally equilibrated Rayleigh–Bénard convection is developed theoretically, with explicit consideration of the role of available potential energy, this being the form in which all the mechanical energy for the flow is supplied. The analysis allows derivation for the first time of a closed analytical expression relating the rate of mixing in symmetric fully developed convection to the rate at which available potential energy is supplied by the thermal forcing. Only about half this supplied energy is dissipated viscously. The remainder is consumed by mixing acting to homogenize the density field. This finding is expected to apply over a wide range of Rayleigh and Prandtl numbers for which the Nusselt number is significantly greater than unity. Thus convection at large Rayleigh number involves energetically efficient mixing of density variations. In contrast to conventional approaches to Rayleigh–Bénard convection, the dissipation of temperature or density variance is shown not to be of direct relevance to the mechanical energy budget. Thus, explicit recognition of available potential energy as the source of mechanical energy for convection, and of both mixing and viscous dissipation as the sinks of this energy, could be of further use in understanding the physics.
Turbulence during the reflection of internal gravity waves at critical and near-critical slopes
- Vamsi K. Chalamalla, Bishakhdatta Gayen, Alberto Scotti, Sutanu Sarkar
-
- Journal:
- Journal of Fluid Mechanics / Volume 729 / 25 August 2013
- Published online by Cambridge University Press:
- 19 July 2013, pp. 47-68
-
- Article
- Export citation
-
Direct numerical simulation is performed with a focus on the characterization of nonlinear dynamics during reflection of a plane internal wave at a sloping bottom. The effect of incoming wave amplitude is assessed by varying the incoming Froude number, $Fr$, and the effect of off-criticality is assessed by varying the slope angle in a range of near-critical values. At low $\mathit{Fr}$, the numerical results agree well with linear inviscid theory of near-critical internal wave reflection. With increasing $\mathit{Fr}$, the reflection process becomes nonlinear with the formation of higher harmonics and the initiation of fine-scale turbulence during the evolution of the reflected wave. Later in time, the wave response becomes quasi-steady with a systematic dependence of turbulence on the temporal and spatial phase. Convective instabilities are found to play a crucial role in the formation of turbulence during each cycle. The cycle evolution of flow statistics is studied in detail and qualitative differences between off-critical and critical reflection are identified. The parametric dependence of turbulence levels on Froude number and slope angle is calculated. Interestingly, at a given value of $\mathit{Fr}$, the turbulent kinetic energy (TKE) can be higher for somewhat off-critical reflection compared to exactly critical reflection. For a fixed slope angle, as the Froude number increases in the simulated cases, the fraction of the input wave energy converted into the turbulent kinetic energy and the fraction of the input wave power dissipated by turbulence also increase.
Horizontal convection dynamics: insights from transient adjustment
- Ross W. Griffiths, Graham O. Hughes, Bishakhdatta Gayen
-
- Journal:
- Journal of Fluid Mechanics / Volume 726 / 10 July 2013
- Published online by Cambridge University Press:
- 11 June 2013, pp. 559-595
-
- Article
- Export citation
-
The dynamics of horizontal convection are revealed by examining transient adjustment toward thermal equilibrium. We restrict attention to high Rayleigh numbers (of $O(1{0}^{12} )$) and a Prandtl number ${\sim }5$ that characterize many practical applications, and consider responses to small changes in the thermal boundary conditions, using laboratory experiments, three-dimensional direct numerical simulations (DNS) and simple theoretical models. The experiments and the mechanical energy budget from the DNS demonstrate that unsteady forcing can produce flow dramatically more active than horizontal convection under steady forcing. The physical mechanisms at work are indicated by the time scales of approach to the new equilibrium, and we show that these can range over two orders of magnitude depending on the imposed change in boundary conditions. Changes that lead to a net destabilizing buoyancy flux give rapid adjustments: for applied heat flux conditions the whole of the circulation is controlled by conduction through the stable portion of the boundary layer, whereas for applied temperature difference the circulation is controlled by small-scale convection within the unstable part of the boundary layer. The experiments, DNS and models are in close agreement and show that the time scale under applied temperatures is as small as 0.01 vertical diffusion time scales, a factor of four smaller than for imposed flux. Both cases give adjustments too rapid for diffusion in the interior to play a significant role, at least through 99 % of the adjustment, and we conclude that diffusion through the full depth is not significant in setting the equilibrium state. Boundary changes leading to a net stabilizing buoyancy flux give a very different response, causing the convection to quickly form a shallow circulation cell, followed eventually by a return to full-depth overturning through a combination of penetrative convection and conduction. The time scale again varies by orders of magnitude, depending on the boundary conditions and the location of the imposed boundary perturbation.
Energetics of horizontal convection
- Bishakhdatta Gayen, Ross W. Griffiths, Graham O. Hughes, Juan A. Saenz
-
- Journal:
- Journal of Fluid Mechanics / Volume 716 / 10 February 2013
- Published online by Cambridge University Press:
- 28 January 2013, R10
-
- Article
- Export citation
-
Three-dimensional direct numerical simulation of horizontal convection is reported for a large Rayleigh number, $\mathit{Ra}\sim O(1{0}^{12} )$, and boundary conditions that allow comparison with previous laboratory experiments. The convection is forced by heating over half of the horizontal base of a long channel and cooling over the other half of the base. The solutions are consistent with the experiments, including small-scale streamwise roll instability developing into a convectively mixed layer within the bottom thermal boundary layer, and a turbulent endwall plume. The mechanical energy budget is shown to be dominated by conversions of available potential energy to kinetic energy by buoyancy flux in the plume and the reverse in the interior of the circulation. These local conversions are three orders of magnitude greater than the total rate of viscous dissipation. The total irreversible mixing is exactly equal to the generation of available potential energy by buoyancy forcing, and one order of magnitude larger than the viscous dissipation. This confirms that dissipation rate is not an indicator of the strength of the circulation and explains why horizontal convection is more energetic than might be expected.
Tidal conversion and turbulence at a model ridge: direct and large eddy simulations
- Narsimha R. Rapaka, Bishakhdatta Gayen, Sutanu Sarkar
-
- Journal:
- Journal of Fluid Mechanics / Volume 715 / 25 January 2013
- Published online by Cambridge University Press:
- 09 January 2013, pp. 181-209
-
- Article
- Export citation
-
Direct and large eddy simulations are performed to study the internal waves generated by the oscillation of a barotropic tide over a model ridge of triangular shape. The objective is to go beyond linear theory and assess the role of nonlinear interactions including turbulence in situations with low tidal excursion number. The criticality parameter, defined as the ratio of the topographic slope to the characteristic slope of the tidal rays, is varied from subcritical to supercritical values. The barotropic tidal forcing is also systematically increased. Numerical results of the energy conversion are compared with linear theory and, in laminar flow at low forcing, they agree well in subcritical and supercritical cases but not at critical slope angle. In critical and supercritical cases with higher forcing, there are convective overturns, turbulence and significant reduction (as much as 25 %) of the radiated wave flux with respect to laminar flow results. Analysis of the baroclinic energy budget and spatial modal analysis are performed to understand the reduction. The near-bottom velocity is intensified at critical angle slope leading to a radiated internal wave beam as well as an upslope bore of cold water with a thermal front. In the critical case, the entire slope has turbulence while, in the supercritical case, turbulence originates near the top of the topography where the slope angle transitions through the critical value. The phase dependence of turbulence within a tidal cycle is examined and found to differ substantially between the ridge slope and the ridge top where the beams from the two sides cross.
Direct and large-eddy simulations of internal tide generation at a near-critical slope
- BISHAKHDATTA GAYEN, SUTANU SARKAR
-
- Journal:
- Journal of Fluid Mechanics / Volume 681 / 25 August 2011
- Published online by Cambridge University Press:
- 25 May 2011, pp. 48-79
-
- Article
- Export citation
-
A numerical study is performed to investigate nonlinear processes during internal wave generation by the oscillation of a background barotropic tide over a sloping bottom. The focus is on the near-critical case where the slope angle is equal to the natural internal wave propagation angle and, consequently, there is a resonant wave response that leads to an intense boundary flow. The resonant wave undergoes both convective and shear instabilities that lead to turbulence with a broad range of scales over the entire slope. A thermal bore is found during upslope flow. Spectra of the baroclinic velocity, both inside the boundary layer and in the external region with free wave propagation, exhibit discrete peaks at the fundamental tidal frequency, higher harmonics of the fundamental, subharmonics and inter-harmonics in addition to a significant continuous part. The internal wave flux and its distribution between the fundamental and harmonics is obtained. Turbulence statistics in the boundary layer including turbulent kinetic energy and dissipation rate are quantified. The slope length is varied with the smaller lengths examined by direct numerical simulation (DNS) and the larger with large-eddy simulation (LES). The peak value of the near-bottom velocity increases with the length of the critical region of the topography. The scaling law that is observed to link the near-bottom peak velocity to slope length is explained by an analytical boundary-layer solution that incorporates an empirically obtained turbulent viscosity. The slope length is also found to have a strong impact on quantities such as the wave energy flux, wave energy spectra, turbulent kinetic energy, turbulent production and turbulent dissipation.
Large eddy simulation of a stratified boundary layer under an oscillatory current
- BISHAKHDATTA GAYEN, SUTANU SARKAR, JOHN R. TAYLOR
-
- Journal:
- Journal of Fluid Mechanics / Volume 643 / 25 January 2010
- Published online by Cambridge University Press:
- 17 December 2009, pp. 233-266
-
- Article
- Export citation
-
A numerical study based on large eddy simulation is performed to investigate a bottom boundary layer under an oscillating tidal current. The focus is on the boundary layer response to an external stratification. The thermal field shows a mixed layer that is separated from the external stratified fluid by a thermocline. The mixed layer grows slowly in time with an oscillatory modulation by the tidal flow. Stratification strongly affects the mean velocity profiles, boundary layer thickness and turbulence levels in the outer region although the effect on the near-bottom unstratified fluid is relatively mild. The turbulence is asymmetric between the accelerating and decelerating stages. The asymmetry is more pronounced with increasing stratification. There is an overshoot of the mean velocity in the outer layer; this jet is linked to the phase asymmetry of the Reynolds shear stress gradient by using the simulation data to examine the mean momentum equation. Depending on the height above the bottom, there is a lag of the maximum turbulent kinetic energy, dissipation and production with respect to the peak external velocity and the value of the lag is found to be influenced by the stratification. Flow instabilities and turbulence in the bottom boundary layer excite internal gravity waves that propagate away into the ambient. Unlike the steady case, the phase lines of the internal waves change direction during the tidal cycle and also from near to far field. The frequency spectrum of the propagating wave field is analysed and found to span a narrow band of frequencies clustered around 45°.
Algebraic and exponential instabilities in a sheared micropolar granular fluid
- BISHAKHDATTA GAYEN, MEHEBOOB
-
- Journal:
- Journal of Fluid Mechanics / Volume 567 / 25 November 2006
- Published online by Cambridge University Press:
- 19 October 2006, pp. 195-233
-
- Article
- Export citation
-
Based on a micropolar continuum of rough granular particles that takes into account the balance equations for the spin (/rotational) velocity and the spin granular temperature, the linear stability characteristics of an unbounded shear flow (${\bm u}\equiv (u_x, u_y, u_z) \,{=}\, (\dot\gamma y, 0, 0)$, where $x$, $y$ and $z$ are the streamwise, transverse and spanwise directions, respectively, and $\dot\gamma$ is the shear rate) are analysed. For pure spanwise perturbations ($k_z\neq 0$, with $k_x\,{=}\,0\,{=}\,k_y$, where $k_i$ is the wavenumber in the $i$th direction), we show that the streamwise translational velocity and the transverse spin velocity modes are subject to linear growths, owing to an inviscid ‘algebraic’ instability (that grows linearly with time). This algebraic instability is shown to be tied to a hidden mechanism of momentum transfer from the translational to the rotational modes, via pure spanwise perturbations to the transverse velocity – in short, we have uncovered an ‘instability-induced rotational-driving’ mechanism. Pure spanwise ($k_z\neq 0$, with $k_x\,{=}\,0\,{=}\,k_y$) and pure transverse ($k_y\neq 0$, with $k_x\,{=}\,0\,{=}\,k_z$) perturbations give rise to ‘exponential’ instabilities (that grow exponentially with time) which are related to similar stationary instabilities in the shear flow of smooth, inelastic particles. Both these instabilities also survive in the limiting case of perfectly elastic but rough particles. The scalings of hydrodynamic modes with wavenumbers have been obtained via the respective long-wave expansion. Perturbations with modulations in all three directions are shown to be stable in the asymptotic time limit, but there could be short-time ‘exponential’ growth of these general perturbations in the long-wave limit for both travelling and stationary waves. The growth rate of all instabilities is maximum at intermediate values of the tangential restitution coefficient ($\beta$), and decreases in both the perfectly smooth ($\beta\to -1$) and rough ($\beta\to 1$) limits; the associated instability length scale is minimum at intermediate $\beta$, and increases in both the perfectly smooth and rough limits. In the perfectly smooth limit, there is a window of particle volume fraction ($\phi$), $\phi_c^s <\phi < \phi_c^t$, over which the flow remains stable to all perturbations. With the inclusion of spin fields, the size of this window decreases and at moderate dissipations with $\beta\,{>}\,0.5$ the flow becomes unstable at all $\phi$.